On the fulfillment of Ward identities in the functional renormalization group approach 
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I consider the fulfillment of conservation laws and Ward 
identities in the one- and two-loop functional renormalization 
group approach. It is shown that in a one-particle irreducible 
scheme of this approach Ward identities are fulfilled only with 
the accuracy of the neglected two-loop terms 0{V^) at one- 
loop order, and with the accuracy 0{V^) at two-loop order 
{Va is the effective interaction vertex at scale A). The one- 
particle self-consistent version of the two-loop RG equations 
which leads to smaller errors in Ward identities due to the 
absence of the terms with non-overlapping loops, is proposed. 
In particular, these modified equations exactly satisfy Ward 
identities in the ladder approximation. 

Strongly-correlated systems attract currently a lot of 
interest. While rigorous mathematical methods, such as 
the Bethe ansatz exist in one dimension, the consider- 
ation of systems with the dimensionality d > 1 requires 
certain approximations. These approximations can break 
conservation laws and therefore lead to results which are 
in contradiction with the conservation of particle num- 
ber, momentum, energy, etc. 

The fulfillment of Ward identities [1-3] guarantees the 
conservation of different quantities. Many years ago a 
way to construct conserving approximations was pro- 
posed by Baym and Kadanoff [4] . They have shown that 
the sufficient condition to fulfill conservation laws is the 
$-derivability, which requires that for a given generating 
functional $[G] the self-energy S and the two-particle 
vertex F irreducible in the particle-hole channel satisfy 
E = S^/SG and T = 6^<i>/6G^. 

Later on, motivated by these ideas, Bickers and 
Scalapino [5] proposed the fluctuation exchange approx- 
imation (FLEX), which is based on a certain choice of 
$ [G] . This approximation has the advantage that it sat- 
isfies the conservation laws by construction. However, 
the disadvantage of FLEX is that the vertex F = STi/SG 
which satisfies Ward identities is obtained only on the 
next level of approximation when the self-energy is al- 
ready computed and therefore the calculation of the self- 
energy is performed with the vertices which do not satisfy 
Ward identities (see e.g. the discussion in Ref. [6]). 

The situation is somewhat better in parquet-type ap- 
proaches [7] (see also Refs. [8,9]). Both, vertex- and self- 
energy effects are accounted for in these approaches on 
the same ground. Although general <i>-derivability of the 
parquet approach was recently proven in Ref. [10], the 
practical application of these approaches for systems with 
d > 1 meets serious computational difficulties and was 
performed only in a restricted number of cases [11,12]. 

Recently proposed functional renormalization group 



(RG) techniques, based on Polchinskii equations [13,14], 
Wick-ordered RG equations [15,16], and one-particle ir- 
reducible (IPI) RG equations [17-19], consider the self- 
energy and vertex corrections on the same footing, but 
require much smaller computational effort. These tech- 
niques were recently applied to the calculation of two- 
particle [14,18-21] and one-particle [22,23] properties of 
the two-dimensional Hubbard model, and to the impurity 
problem in a Luttinger liquid [24] . 

The fulfillment of Ward identities within the functional 
RG techniques is a non-trivial problem because of the use 
of different approximations: the truncation of an infinite 
hierarchy of RG equations (loop expansion), projection of 
the vertices on the Fermi surface etc. The main restric- 
tion comes from the loop expansion, used within these 
approaches. Unlike previous field-theoretical and Wilson 
RG approaches, not only singular, but also regular terms 
are accounted for in functional RG. This can lead to some 
inconsistencies and to the violation of Ward identities. 
Besides the approximations mentioned above, the intro- 
duction of the cutoffs in momentum or energy spaces can 
break Ward identities at the intermediate stages of the 
RG flow [25]. 

Therefore, it is useful to verify whether Ward identi- 
ties are fulflUed in the functional RG approach. In the 
present paper I consider the fulfillment of Ward identities 
in one- and two-loop IPI functional RG analysis without 
the projection of the vertices to the Fermi surface. 

I consider a model with the action 
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with the general (frequency dependent) interaction 
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The conservation of charge leads to a continuity equa- 
tion, which can be written in the continuum limit as 



q^J^ = 0, 
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where /i = 0...3, q — {u!,q), Jq* is the contribution of 
electron-electron interaction to the current, Jq' = for 
a density-density interaction with y^^^j,^.^ j.^ — f{ki — k^). 
The Ward identity which follows from (2) connects the 
one-particle self-energy with the charge- and charge- 
current vertices. 
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where index ir stands for connected contributions, Gfe is 
full one-particle Green function, r^j^^. is the contribution 
to charge-current vertex produced by Jq*. The corre- 
sponding Ward identity for the model (1) reads (see Refs. 
[2,3,26] for a derivation) 

y^(~^^fc.fc':fc-l-n.fc'-n + ^k,k' 'M' -q,k+q) 
k' 

xGk'Gk'-q = Sfe - Sfe+g, (4) 

where = (ro,r). 

The functional renormalization-group approach leads 
to an infinite hierarchy of differential equations for the 
self-energy E, two-particle vertex V^'^^ = V, and higher- 
order vertices V^"'\ n > 4. Truncated at the one-loop 
order, i.e. after neglect of y^"^ with n > 4 this hierar- 
chy reduces to two intcgro-differential equations for the 
self-energy and the two-particle vertex, which in IPI RG 
scheme [17,18] have the form 



'dA 
'dA 



= o S^, (5a) 



o (G^ oS^ + S^o G^) o V^. (5b) 



In Eq. (5) o stands for the summation over intermedi- 
ate momentum-, frequency and spin- variables according 
to standard diagrammatic rules, A is the scaling param- 
eter specified below. The initial condition for Eqs. (5) is 
T.^" = 0, V^" = Vo where Aq > max(|£fe]). The propa- 
gators S, G are connected by 

5^ = -(G'^r^iG^)-', 

where Gg is the non-interacting Green function at scale 

A. In particular, in the sharp momentum cutoff scheme 
[18,22] the propagators G^ and are given by 



Gt , 

St 



f?(kk|-A) 

-'^(lekl - A) / ii.„-ek 



(6) 



in the sharp frequency cutoff scheme (at T = 0) [24] by 



Gt \ 

St - 



(7) 



and in the temperature cutoff scheme [19] A — T and 
T'/' .T_ 1 (»^n + £k)/TV^ 



Note that the results obtained in the present paper do 
not depend on the sharpness of the cutoff procedure in 



momentum or energy space and are valid for smooth cut- 
off functions as well. 

Consider first the Ward identity (4) for pure frequency 
shift q = (w, 0) within the momentum-cutoff RG ap- 
proach (6). Written in differential form the correspond- 
ing identity reads 



dT,k/dko = '^{Vk,k';k',k - '^Vk,k';k,k')Gl, 



(9) 



where k = {ko, k). 

To verify whether this identity is satisfied at each in- 
termediate scale A within the 1-loop RG scheme (5), one 
can use the following procedure: differentiate (9) with 
respect to A and compare with RG equations (5) differ- 
entiated with respect to ko- The following identities are 
helpful when performing this procedure: 
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In this way one finds that the identity (9) is satisfied if 
and only if 

dVk 

~ ^ 2_^p ^k,k-\-p;k-\-p,k^k-\-p,k+q;k+q,k-\-p^k-\-p 

Vk ,k-\-q-\-p;k-\-q,k-\-p ) (11) 
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+ yk,k+p:k,k+p{yk+p.k+q;k+q,k+p — '^yk+p.k+q;k+p,k+q)\ 
yk,k+p+q;k+p,k+q{Gk+p+qGk+p + Gk+p+qGk+p) 



+ yk,k+q;k+p,q-p+k {Gk+pGq-p+k + Gg_p_^_kGk+p) 

Note that these equations have the same structure as the 
RG equation for the vertex, Eq. (5b) with the replace- 
ments dV/dA dV/dko and S ^ G'^, and therefore can 
be shortly written as 



dko 



o [G^ o (G^)^ + (G^)^ o Ga] o Va (12) 



To see whether this equation is compatible with the RG 
equation for the vertex, Eq. (5b), one can differentiate 
Eq. (11) once more with respect to A and compare it 
with (5b) differentiated with respect to ko- After long 
algebraic manipulations one obtains for the difference 
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Eq. (5b) 
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(and similar for VjAfc+^.j, j^^^) two types of terms, which 
can be written in a short notation as 



So 



(13) 



Both terms in Eq. (13) are of the order 0{V^) and closer 
inspection of these terms shows that the terms of each 
type do not cancel each other. Therefore, the momentum 
cutoff RG approach in one-loop approximation does not 
satisfy the Ward identity (9) exactly, but rather with the 
accuracy O(V^). This accuracy is the same as the accu- 
racy of the one-loop equations (5). The terms in Eq. (13) 
have a different structure: in the terminology of Ref. [15] 
the terms in the first and second lines of Eq. (13) corre- 
spond to diagrams with non-overlapping and overlapping 
loops, respectively. The presence of the terms with non- 
overlapping loops prevents the use of more strict power 
counting arguments [15] to find the bounds for the cor- 
responding contributions responsible for the violation of 
the Ward identity. 





a) 



b) 



FIG. 1. The two-loop corrections to the two-particle ver- 
tex: (a) vertex correction, (b) self-energy correction. The 
lines with one and two cutting slashes correspond to the 
and S"^ propagators, other lines to G"^ propagators. 

Consider now the two-loop IPI RG approach. The 
two-loop approximation corresponds to the inclusion of 
the contribution of the three-particle vertex V^^^ up to 
the order O(V^). The corresponding equations read 



V''o{G^oS'' + S^oG^)oV^ 



(14a) 
(14b) 
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J A 

where the last two lines of Eq. (14b) correspond to the 
vertex (Fig. la) and self-energy (Fig. lb) corrections 
to the one-loop equations. Performing the same steps as 
for the one-loop approximation, we obtain a necessary 
condition to fulfill Ward identities: 



OKq OKo 



(15) 
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(G^)^ O / dA'iV"^ O G^ O ^qA ^yA ^gA^ 

A 
An 
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^(G^y o / dt^{V^ o {G'' Y o o o ) 

J A 

Comparing again this equation with the RG equation for 
the vertex, Eq. (16b), one can observe that the equations 
are compatible with the accuracy O(V^). As well as in 
the one-loop approximation, both types of terms, with 
overlapping and non-overlapping loops of order O(V^), 
violate the Ward identity (9). 

However the modification of Eqs. (14) with the last 
term replaced by its local-in-A analog, i.e. 



S'^ o / dA'iV^ o G^ oV^ o G^ oV^ oS^) 

I A 

S^oV^oG^oV^oG^oV^o G^ 



leads to improvement in the fuUfilment of the Ward iden- 
tities. The corresponding equations read 
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(16b) 
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The difference between (14) and (16) is itself of the or- 
der of 0(F*), i.e. of the order of the terms neglected in 
the two-loop approximation. The advantage of the modi- 
fied two-loop RG equations, Eq. (16) is that they violate 
the Ward identity (9) by terms with overlapping loops 
only. Since such terms can correct only vertices but not 
the internal Green functions, one can consider Eqs. (16) 
as one-particle self-consistent. The absence of the terms 
with non-overlapping loops in the difference between the 
r.h.s. and the l.h.s. of the Ward identity (9) in this case 
is a consequence of the non-trivial cancellation between 
derivatives of the self-energy contained in the propaga- 
tors of IPI scheme (6) and the vertex corrections. 

In particular, in the ladder approximation, when one 
selects only diagrams in the appropriate particle-hole 
channel and neglects the last term in the r.h.s. of Eq. 
(16b), Eqs. (16) take the form 



"dA 

dVA 
dA 



(17a) 



^ = t^'' ° (G"" ° ° ^''Iph, ladder (17b) 

dA 



These ladder equations differ from those in the one- 
loop approximation (5) only by the replacement S'a — > 
dG^/dA in the equation for the two-particle vertex and 
fulfill Ward identities exactly at any scale A. The solu- 
tion of these equations leads to RPA-type vertices Va, 
which satisfy standard integral equations 
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kp;pk kp;pk ' / ^ * kk' ]k' k^y k' ) ^k'p\pk' 

k' 

=v° -l\^v° (r^ )'^v^ 
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+ ^^i^kk';k'k ~ "^^kk' M'^ ^^k')"^^!^ p;k' p (18) 
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and the mean-ficld-typc result for the self-energy, 

^k = y^(^fcp;pfc ~ '^^kp;kp)^p 



(19) 



with the full (dressed) Green function G^. In the limit 
A — > Eqs. (18) and (19) reduce to standard RPA and 
mean-field results, respectively. 

For the non-ladder case and the regular curved Fermi 
surface one can use power counting arguments [15,17] to 
show that the equations (16) violate the Ward identity 
(9) by finite (non-divergent in the A ^ limit) terms 
O(V^). One can expect that these equations will bet- 
ter satisfy Ward identities for both regular and singular 
Fermi surfaces due to the smallness of the contributions 
of the terms with overlapping loops [15]. Note that the 
replacement of the last term in Eq. (16) by its local ana- 
log and the additional replacement Sa dGA/dA in the 
same term (which again produces finite errors of the or- 
der 0{V^)) can further improve the fulfilllment of Ward 
identities. The analysis of this possibility requires how- 
ever the consideration of the three-loop terms, which is 
beyond the scope of the present paper. 

Therefore, in general, the momentum-cutoff RG 
scheme (6) at two-loop order satisfies the frequency com- 
ponent of the Ward identity (4) with the accuracy O(V^). 
The same situation holds for the momentum component 
of the Ward identity (4) in the frequency cutoff RG ap- 
proach, Eq. (7). The frequency component of the Ward 
identity (4) in frequency cutoff scheme and momentum 
component in momentum-cutoff scheme are violated be- 
cause of additional terms coming from the differentiation 
of the cutoff functions [25]. At the same time, in the 
temperature cutoff RG scheme (8), which does not con- 
tain any cutoff in momentum or frequency space, similar 
calculations show that both, frequency- and momentum- 
components of the Ward identity are satisfied at two-loop 
order up to the same accuracy 0{V^). 

In conclusion, we have considered the fulfillment of 
Ward identities in the one- and two-loop RG approach. 
The frequency components of Ward identities in the 
momentum-cutoff scheme and the momentum compo- 
nents in the frequency-cutoff scheme are generally sat- 
isfied up to terms of O(V^) in the one-loop approxima- 
tion and O(V^) in the two-loop approximation. The 
one-particle self-consistent version of the two-loop RG 
equations (16) leads to smaller errors in the fulfillment 
of Ward identities due to absence of terms with non- 
overlapping loops. In the ladder approximation this ver- 
sion of the two-loop equations satisfies Ward identities 
exactly. 



[1] J. C. Ward, Phys. Rev. 78, 182 (1950). 

[2] A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinskii, 

Methods of quantum field theory in statistical physics, 

Prentice-Hall, New-Jercey, 1963, §19. 
[3] T. Toyoda, Ann. Phys. 173, 226 (1997) 
[4] G. Baym and L. P. Kadanoff, Phys. Rev. 124, 287 

(1961). 

[5] N. E. Bickers, D. J. Scalapino, and S. R. White, Phys. 

Rev. Lett. 62, 961 (1989); N. E. Bickers and D. J. 

Scalapino, Ann. Phys. 193, 206 (1989). 
[6] Y. M. Vilk and A. M. S. Tremblay, J. Phys. I (Prance) 

7, 1309 (1997). 

[7] A. A. Abrikosov, Physics 2, 5 (1965); Yu. A. Bychkov, 
L. P. Gor'kov, and I. E. Dzyaloshinskii, Zh. Eksp. Teor. 
Fiz. 50, 738 (1966) [Sov. Phys. JETP, 23, 489 (1966)] 

[8] N. E. Bickers and S. R. White, Phys. Rev. B 43, 8044 
(1991). 

[9] V. Janis, J. Phys.: Cond. Matt. 10, 2915 (1998). 
[10] V. Janis, cond-mat/9806118 (unpublished). 
[11] V. Janis, Phys. Rev. B 60, 11345 (1999). 
[12] I. E. Dzyaloshinskii and V. M. Yakovcnko, Zh. Eksp. 

Teor. Fiz. 94, 344 (1988) [Sov. Phys. JETP 67, 844 

(1988)]; Int. J. Mod. Phys. B 2, 667 (1988); A. T. 

Zholcznyak, V. M. Yakovenko, and I. E. Dzyaloshinskii, 

Phys. Rev. B 55, 3200 (1997). 
[13] J. Polchinski, Nucl. Phys. B 231, 269 (1984). 
[14] D. Zanchi and H. J. Schulz, Phys. Rev. B 54, 9509 (1996); 

ibid. 61, 13609 (2000). 
[15] M. Salmhofcr, Rcnormalization, Springer, Berlin, 1999. 
[16] C. J. Halboth and W. Metzner, Phys. Rev. Lett. 85, 5162 

(2000) ; Phys. Rev. B 61, 7364 (2000). 

[17] M. Salmhofer and C. Honerkamp, Progr. Theor. Phys. 
105, 1 (2001). 

[18] C. Honerkamp, M. Salmhofcr, N. Furukawa, and T. M. 

Rice, Phys. Rev. B 63, 035109 (2001). 
[19] C. Honerkamp and M. Salmhofer, Phys. Rev. Lett. 87 

(2001) 187004; Phys. Rev. B 64, 184516 (2001). 
[20] C. Honerkamp, Eur. Phys. J. B 21, 81 (2001). 

[21] A. P. Kampf and A. A. Katanin, Phys. Rev. B 67, 125104 
(2003); A. A. Katanin and A. P. Kampf, Phys. Rev. B 
68, 195101 (2003). 

[22] C. Honerkamp and M. Salmhofer, Phys. Rev. B 67, 
174504 (2003). 

[23] A. A. Katanin and A. P. Kampf, cond-mat/0310112 (un- 
published). 

[24] V. Meden, W. Metzner, U. Schollwock, and K. 

Schonhammor, Phys. Rev. B 65, 045318 (2002). 
[25] T. Enss, PhD thesis, in preparation (unpublished). 
[26] M. Rcvzcn, T. Toyoda, Y. Takahashi, and F. C. Khanna, 

Phys. Rev. B 40, 769771 (1989). 



4 



